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Abstract 

We study the propagation of ultra-short short sohtons in a cubic nonhnear medium modeled 
by nonlinear Maxwell's equations with stochastic variations of media. We consider three cases: 
variations of (a) the dispersion, (b) the phase velocity, (c) the nonlinear coefficient. Using a mod- 
ified multi-scale expansion for stochastic systems, we derive new stochastic generalizations of the 
short pulse equation that approximate the solutions of stochastic nonlinear Maxwell's equations. 
Numerical simulations show that soliton solutions of the short pulse equation propagate stably in 
stochastic nonlinear Maxwell's equations and that the generalized stochastic short pulse equations 
approximate the solutions to the stochastic Maxwell's equations over the distances under consid- 
eration. This holds for both a pathwise comparison of the stochastic equations as well as for a 
comparison of the resulting probability densities. 
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I. INTRODUCTION 



Many integrable partial differential equations arise as asymptotic expansions of nonlinear 
wave equations. In the particular context of ultra-short pulses propagating in cubic nonlinear 
media, the basic model equation is given by Maxwell's wave equation with dispersion and 
nonlinearity coming from the material's response to the excitation by the electric field [1]. 
When the pulse width is large in comparison to the carrier wavelength, an asymptotic 
expansion of the nonlinear wave equation yields the cubic nonlinear Schrddinger equation 
(NLSE) [2]. The NLSE is integrable [3] and has bright soliton solutions that have been 
studied intensively over the last four decades, not only because of their mathematical beauty, 
but also because of a variety of possible applications in modern optical technology P| . Recent 
experimental progress in the field of ultra-short pulses [3] has led to increased interest in 
models that go beyond the NLSE description: As the length of the pulse shortens, the basic 
assumption of scale separation between envelope and carrier wave is not satisfied anymore 
[6] and additional correction terms have to be added [7j. On the other hand, it is possible 
to introduce a different scaling [8j in the asymptotic analysis of the nonlinear wave equation 
that is appropriate for ultra-short pulses. In this way, it is possible to derive a different 
approximation of the solution of Maxwell's equation, the short pulse equation (SPE) [H]. 
Although derived under entirely different assumptions in comparison to the NLSE, the SPE 
is integrable [10] as well and was shown to possess bright solitons that may be as short as 
three cycles [H]. Consequently, over the last years, the mathematical structure of the SPE 
was subject to intensive research [T2HT7j. 

When a solitary wave solution is taken as the initial condition in the SPE, it persists for 
arbitrary propagation distances. As initial condition of the original nonlinear wave equation, 
however, it is expected to change its shape due to the fact that the SPE is only an approx- 
imation to the nonlinear wave equation. If e is the expansion factor in the multiple-scale 
asymptotic analysis, we would expect to see a noticeable deviation only after propagation 
distances of the order of 1/e [TS]. We show at the end of the following section numerically 
by propagating the soliton to a distance of order that this change happens indeed very 
slowly, such that we can say that the SPE solitons persist in Maxwell's equations. The 
slow change of the soliton solution is governed by higher-order effects that appear as addi- 
tional terms in the multi-scale expansion and that can be studied systematically [19]. In the 
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present paper, since we are focusing on the leading order impact of stochastic perturbations 
of the parameters in the equations, we neglect these higher-order deterministic effects. In- 
corporating deterministic higher-order effects in the presence of random perturbations lies 
beyond the scope of this paper and will be studied in the future. 

Such random fluctuations are widely present in nature and models in which stochasticity 
is taken into account are more realistic than deterministic models. One important question 
that arises in the context of integrable systems is whether soliton solutions persist when 
perturbed stochastically: NLSE solitons broaden under the influence of stochastic variations 
of the dispersion coefficient [20] . For the NLSE case, however, this broadening is sufficiently 
slow such that the soliton still persists over long distances, provided that the stochastic 
perturbations are sufficiently small. One result of this paper is that SPE solitons show 
similar behavior. 

When studying the impact of stochastic variations of the media, one faces the problem of 
coarse- graining noise, namely the question of how microscopic variations of the parameters in 
the nonlinear wave equation relate to variations of the parameters in the SPE or the NLSE. 



In section III , based on previous work [211 122] , we show how this coarse-graining of noise can 
be done explicitly in the multi-scale expansion that leads from the nonlinear wave equation 
to the SPE: We introduce randomness in the linear dispersive part of the susceptibility of 
the media and obtain a stochastic dispersion coefficient in the SPE, yielding a stochastic 
generalization of the short pulse equation. Similarly, stochasticity can be introduced in the 
phase velocity and and in the nonlinearity. The resulting generalized stochastic short pulse 
equations constitute the main result of this work. 

It remains to be shown that the solutions to the stochastic short pulse equations stay close 
to the solutions of the original stochastic nonlinear wave equation. Even in the deterministic 
context this is very challenging and, in fact, this problem has only be solved recently |23] . 



Therefore, in section IV we present a numerical comparison of the solutions to the nonlinear 
wave equation with randomness and the corresponding stochastic short pulse equations. We 
show that, over the distances under consideration, SPE solitons also persist in the stochastic 
short pulse equations and in the stochastic nonlinear wave equation. Moreover, extensive 
numerical simulations show that the SPE model can be used to predict the likelihood of large 
deviations from the deterministic solution when we compare particular stochastic paths and 
as well as the corresponding probability distributions. 



II. ULTRA-SHORT SOLITONS IN DETERMINISTIC MAXWELL'S EQUATIONS 

In this section we briefly review the setting of the short pulse equation and its derivation 
from Maxwell's equations. We flrst consider the following deterministic one-dimensional 
wave equation describing pulse propagation in bulk silica [9] 



Here, x is our evolution variable and we assume that we are given initial conditions u{x = 0,t) 
and Ux{x = 0,t). In this section we assume the coefficients xo and xs to be fixed, in the 
next section we will modify the problem by allowing the coefficients xo and xs to vary 
stochastically. Above equation can be directly derived from Maxwell's equations under 
the assumption that the response of the medium is instantaneous, otherwise a retarded 
response of the material will lead to convolution integrals in the polarization which will 
yield a nonlocal wave equation [23] . We also assume that the Fourier transform u is zero for 
very low wavelengths. This accounts for the fact that such waves cannot propagate in the 
optical medium. 

As for many wave equations, we can use asymptotic expansions in order to construct 
approximate solutions to ([T|. The short pulse equation (SPE) can be derived from ([T]) 
using a multiple-scale expansion |i9j or an equivalent method like the renormalization group 
method [2S]: A multi-scale expansion of the form 




(1) 




(2) 



with 



t — X 



(3) 



e 



yields as leading order 



2d^dx,Ao = XoAo + X3dlAl . 



(4) 



This form of the short pulse equation can be transformed to the standard form 




XX 



(5) 



through the transformation given by 




X 




(6) 
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Note that we change here from the physical variables (xi, 0) (which are rescaled and shifted 
with respect to the variables {x,t) of the nonlinear wave equation) to the mathematical 
variables (T,X). In all simulation results later, we will always use the physical variables 
(x, t) of the original nonlinear wave equation. 

To our knowledge, it is unknown whether the original wave equation ([T| is integrable. 
Sakovich and Sakovich, however, have shown that the SPE (|5| is integrable and, moreover, 
admits one-soliton solutions [11] that can be written as 

m sin ip sinh + n cos ip cosh 



U = Amn- 



m? sin^ -0 + cosh^ 



m sin 2'?/' — n sinh 20 

X = Y + 2mn ^ . , tt^ (7) 

sm ip + n'^ cosh 

with 

= m(F + T), i) = n{Y-T), n = ^/il-m"^). (8) 

This family of solutions depends on the parameter m, and the condition for these solitary 
wave solutions being non-singular was found to be 

m < rricr = sin - ^ 0.383. (9) 
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The parameter m determines the shortness of the pulse. As m reaches its critical value, the 
shortest pulse of this family we can construct is approximately a three cycle pulse. When 
m is small (for instance, m = 0.01), solitary wave solution of the SPE can be approximated 

by 

t/^4m cos(X-T)sech(m(X + T)) (10) 

such that the pulse is similar to a bright soliton solution of the cubic nonlinear Schrodinger 
equation. In the present work, we focus on the behavior of ultra-short solutions. In our 
simulations, we take m = 0.3. 

With these soliton solutions of the SPE at hand, the question arises whether they also 
show stable propagation when taken as initial conditions in the original problem ([T]). First, 
it was shown for the linear case (xa = 0) that the solution of the SPE stays close to the 
solution of the original problem [25]. For the nonlinear case, the resulting proof is much 
more difficult and a solution was found only recently [23] . 

Figure [l] shows a numerical solution of the nonlinear wave equation ([T| with the SPE 
soliton solution taken as initial conditions. In the simulation, we chose the parameter m 



to be 0.3 and e = 0.2. This snapshot of the sohton profile is compared to the analytical 
solution of the SPE given by ([T]). The soliton shows stable propagation in the nonlinear wave 
equation - note that the evolution distance in x is 25.625 ~ (9(l/e^), which is significantly 
further than we would expect for a leading order multi-scale expansion. The reason for this 
can be found when considering addtional higher-order terms in the multi-scale expansion 
[19] . In this work, however, we neglect such high-order contributions. 




FIG. 1. Profile of the SPE soliton in the nonlinear wave equation (|T]). The solid line shows 
the snapshot of the soliton solution as a result of its propagation in Maxwell's nonlinear wave 
equation. The dashed line is the corresponding solution to the short pulse equation with the same 
initial condition. The small deviations are due to higher-order deterministic corrections. 

In general, we expect the difference between the solution u of the nonlinear wave equation 
and the soliton Us to grow with the propagation distance x. We can quantify this growth 
by a numerical comparison. In figure |2] we plot the evolution of several norms of u — Ug, in 
particular the L^, L"^ and L°°-norms defined by 

\\f\\L^=(^J \fittdty\ ||/|U.=max(|/(t)|). (11) 

Whereas the and the L^-norms display monotonic growth, the L°°-norm shows oscilla- 
tions that are due to the fact that this norm is very sensitive to shifts in time as it compares 
the solutions pointwise. 
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FIG. 2. Evolution of the deviations of the SPE sohton when propagating in the nonhnear wave 
equation ([l]). The sohd line shows the evolution of the L^-norm given by \\u — Us\\li, the dashed 
line the evolution of the L°^-norm, and the dotted line the evolution of the L^-norm. 

III. THE STOCHASTIC SHORT PULSE EQUATION 

Optical soliton propagation in fibers in the presence of a stochastic perturbation has been 
studied intensively in the context of the NLSE model [20]. The sources of randomness in 
optical fibers may vary. For example, the phase of the wave may fluctuate due to the presence 
of the randomness. One major source of medium-related stochastic phase fluctuations are 
the inhomogeneities in a flber's core or the fluctuations in the linear refractive index of the 
core. Stimulated Brillouin scattering, stimulated Raman scattering, and inhomogeneities of 
the medium are the possible sources of the phase fluctuations [261 EZl- The stochasticity 
may come from the nonlinearity of the medium as well. The dynamical effect of the noise 
added by the stochastic nature of the nonlinearity is much smaller than the noise due to 
the inhomogeneities [2B] . Nevertheless, there may be other sources of randomness playing a 
role, and they may originate from the other parts of the system such as the inherent power 
fluctuation in lasers used as input pumps. Apart from these, quantum phase fluctuations 
are also sources of phase noise in optical flbers although they are practically negligible [7]. 
In many applications, Langevin noise is used to study the fluctuations in the system [1]. 
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In the present work, we do not model the precise microscopic origins of the randomness, 
but rather assume that they lead to small stochastic variations of the coefficients of the 
equations. This approach is common, one particular case studied in the past is the NLSE 
with a linear multiplicative stochastic term corresponding to stochastic variations of the 
dispersion. In analogy, it is natural to ask how the ultra-short solitons of the deterministic 
SPE propagate in a stochastic environment. 

Therefore, the simplest corresponding stochastic variations can be introduced to the wave 
equation ([T]) by making the coefficient Xo representing the linear susceptibility stochastic. 
Setting 

XO = XO + Xrand(a;), Xrand(a;) = l^^{x) , (12) 

leads to the stochastic wave equation such that 

Uxx - utt = [Xo + t^^ix)] u + x?.{u^)tt , (13) 

where represents white noise, hence (^(a;)^(x')) = 5{x — x'), and the coefficient u measures 
the strength of the noise. Although mathematically convenient, we note that this is only an 
approximation, as in reality, the noise will be short-correlated rather than delta-correlated. 
In the numerical simulations we also will assume that the coefficient of the dispersive term 
Xo will always be larger than zero which correctly represents the underlying physical system. 

In order to derive a stochastic generalization of the short pulse equation, we introduce 
the same multi-scale expansion as in the deterministic case of the form 

u{x,t) = eMo(0,xo,xi,X2, ...) + e^Mi(0, xq, xi, 0:2, ...) + ... (14) 

with again rescaled variables 

(j) = , Xn = (15) 

As in the deterministic case, we solve now the equation order by order. As terms of the 
order 1/e cancel, we obtain at 0{1) the equation Mq^^^ = such that the leading order Mq 
is independent of the fast variable Xq. The equation for Mi is given by a combination of 
terms and can be written as: 

- 2(Mi),„^ = 2(Mo),,^ + (xo + '^^{x))Mo + XsiMl)^ ■ (16) 
In the deterministic case, if u = 0, the solvability condition of this equation yields the short 



pulse equation, and therefore the entire r.h.s. of (16) will be zero. For the stochastic case, 



the standard multi-scale argument has to be slightly modified due to the presence of the 



noise term. If we integrate (16) with respect to the fast variable Xq from xq = to a;o = 1, 
the noise term ^(x) yields a random number (q given by Co = Jq ^{x)dx. As ^ is white noise, 
Co is normally distributed with mean zero and variance one. On the slow scale xi, the points 
and 1 correspond to and e and Co represents the cumulative microscopic noise for this 
interval. For the next interval [1, 2) on the fast scale (corresponding to [e, 2e) on the slow 
scale) we can repeat the process and obtain d = ^{x)dx. In this way we obtain a sequence 
of independent normally distributed random numbers {(k) = (Co;Ci; •••)• ^^e slow scale, 
the random number characterizes the noise in the fc-th interval [ke, {k + l)e). Hence, we 
can view the sequence of random numbers {(k) on the slow scale as a discretization of a 
white noise process H = with {E{xi)E{x[)) = e6{xi — x[). In this way we obtain the 

stochastic short pulse equation with coarse-grained noise as 

- 2(Mo).,0 = (xo + i^H(xi)) Mo + X3{M^)4>^ (17) 
One important difference between the deterministic and the stochastic case is that, if we 



introduce (17) back in (16), we obtain a dependence of Mi on the fast variable Xq that is not 
present in the deterministic case. This is especially important when considering the next 
order 0{e^) in the asymptotic expansion. In this work, however, we only consider terms up 
to the order 0{e). 

Note that we obtain the strength of the slow-scale noise S explicitly and that both the 
microscopic strength of the noise (given by u) and the expansion coefficient e infiuence the 
strength of the coarse-grained noise. The reason for the fact that the coarse-graining of the 
noise takes a rather simple form lies in the fact that the randomness is introduced in the 
linear dispersive term and that the leading order evolution is trivial as Mo is independent of 
Xo- In cases where the leading order solution is more complicated, its structure will influence 
the strength of the coarse-grained noise [22] . 

Aside from fluctuations in the dispersion, other parameters might vary stochastically as 
well. In deriving the short pulse equation from Maxwell's equations, an approximate flt of 
the linear susceptibility of the form 

X^'K\)-^xi'^-X?\' (18) 

in Fourier domain [9l |29] was made. This susceptibility approximation is found by fltting 
the experimental data for light propagation in silica flbers in the infrared regime. The flrst 



stochastic generalization of the SPE treated above is derived by considering the fluctuations 
in the A-dependent term. However, one can argue that the fluctuations may come from the 
rather than the X2^'* term, corresponding to stochastic variations of the 
phase velocity. This treatment will lead to a different stochastic wave equation: 

Uxx - [1 + ^^(x)] Utt = Xou + X3iu^)tt ■ (19) 

As before, a multiple scale expansion, see ^ and ([s]), leads to a different form of the 
stochastic SPE. As the calculations are similar to the previous case of random dispersion, 
we just give the result: 

- 2(Mo).,<^ = XoMo + X3{M^)h + (^S(^i)) (^0)00 (20) 

Finally, we will introduce the noise in the nonlinearity as well, which will result in a 
stochastic wave equation in the form of 

Uxx - Utt = XoU + [X3 + t^^i^)] {U^)tt ■ (21) 

The corresponding stochastic SPE takes the form 

- 2(Mo),,^ = XoMo + (X3 + '^^{xi)) {MD^^ . (22) 

If one sets z/ = in any of these stochastic generalizations of the SPE, one will obtain the 
deterministic short pulse equation (|4]). 



X"q^ term in (18 



IV. NUMERICAL COMPARISON 



A. Pathwise comparison 

Again, we expect (17) to be a decent approximation of the nonlinear wave equation ([T]) 



with a stochastic linear susceptibility given by (12) for propagation distances that are of the 
order of at least 0{l/e). Since we obtained a very good agreement up to distances of C(l/e^) 
in the deterministic case, we will consider this as our maximum propagation distance for the 
stochastic well. 

In order to evaluate the quality of the approximation numerically, we can look at a path- 
wise comparison. From the coarse-graining procedure described above, we obtain for any 
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particular realization of the micro-scale noise ^ on the fast scale a macroscopic noise Sp 



on the slow scale. Then, we can compare the solution of the stochastic SPE (17) with the 



solution of the nonlinear wave equation for this particular realization of the noise. Again we 
choose to look at the differences of the solution of the stochastic PDEs and the deterministic 
soliton. For the simulation of the nonlinear wave equation, we subtract the error that exists 
in the deterministic system, presented in figure [2j As a result, we obtain two stochastic 
processes - one given by the difference of the solution of the stochastic SPE and the de- 
terministic soliton and the other given by the difference of the solution of the stochastic 
nonlinear wave equation and the deterministic solution to the wave equation (with the SPE 
soliton as initial condition). 

The result of the corresponding numerical simulations is shown in figure |3] and shows 
excellent agreement of the approximation and the solution of the original equation. 
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FIG. 3. Stochastic dispersion: Evolution of the stochastic difference of the SPE soHton and the 
solution of the nonlinear wave equation ([T]). The solid line shows the evolution of the L°°-norm of 
the difference generated by the solution of the stochastic nonlinear wave equation and the dashed 
line shows the corresponding difference for the corresponding evolution in the stochastic short pulse 
equation. The strength parameter of the microscopic noise was 0.125. 



The next figure |4] shows the same comparison for the case of randomness in the phase 
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velocity. In this case, we compare the solutions to (20) to the solutions of ([T]) for the case 
of random variations of Xo^^- 
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FIG. 4. Stochastic phase velocity: Evolution of the stochastic difference of the SPE soliton and the 
solution of the nonlinear wave equation ([T]) for the case of random phase velocity. The solid line 
shows the evolution of the L°°-norm of the difference generated by the solution of the stochastic 
nonlinear wave equation and the dashed line shows the corresponding difference for the correspond- 
ing evolution in the stochastic short pulse equation. The strength parameter v of the microscopic 
noise was 0.005. Recall that e = 0.2 such that the ratio v/e^ is still fairly small. 



The third case is the comparison for the case of random fluctuation of the nonlinear 
coefficient. Figure [5] compares the solution of the stochastic short pulse equation with the 
solution to stochastic Maxwell's equations. 

Hence, for all three cases, we have numerical evidence that the coarse-grained stochastic 
generalizations of the SPE capture correctly the impact of the noise on the solitons when 
propagating in the stochastic nonlinear wave equation. 
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FIG. 5. Stochastic phase velocity: Evolution of the stochastic difference of the SPE soliton and the 
solution of the nonlinear wave equation ([T]) for the case of random phase velocity. The solid line 
shows the evolution of the L°°-norm of the difference generated by the solution of the stochastic 
nonlinear wave equation and the dashed line shows the corresponding difference for the correspond- 
ing evolution in the stochastic short pulse equation. The strength parameter v of the microscopic 
noise was 0.05. 

B. Probability distribution 

This pathwise correspondence indicates that we can also expect a decent agreement in the 
prediction of the related probability distributions. To check this, we can plot the histograms 
of the deviations of the solutions to the stochastic PDEs from the deterministic evolution of 
the SPE solitons. The following figure [6] was obtained by creating a probability distributions 
from 10000 realizations of the stochastic short pulse equation and stochastic Maxwell's 
equation for the case of stochastic variations of the dispersion. 

Although we see a small difference of the SPE and stochastic Maxwell's equations for small 
deviations, we obtain overall an excellent agreement between the approximation given by the 
stochastic short pulse equation and the stochastic nonlinear wave equation. In particular, 
the shape of the probability distribution is correctly predicted by the SPE approximation 
and large deviations from the mean are correctly accounted for as well. 
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FIG. 6. Probability distributions of the stochastic difference of the SPE sohton and the solution 
of the nonlinear wave equation ([T]) at the propagation distance x = 25.625. The x-axis shows 
the L°°-norm of the difference generated by the solution of the stochastic PDEs in comparison 
to the deterministic solution (solid line: Maxwell, dashes: SPE). In this way, the solution of the 
nonlinear wave equation has been corrected for the deterministic error. The strength of the noise 
was = 0.05. 



V. CONCLUSION 



In this paper we have derived stochastic generalizations of the short pulse equation and 
shown explicitly how the noise processes governing the parameters of the SPE can be con- 
structed from the noise present in the underlying nonlinear wave equation. Moreover, we 
have shown that SPE solitons persist in the nonlinear wave equation with stochastic vari- 
ations present and that they exhibit only small changes over the propagation distances 
under consideration. The changes arising from the random part of the susceptibility can be 
correctly characterized by the stochastic short pulse equation. 
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